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Abstract. A catalogue of light curve solutions of contact binary stars has been compiled. It contains the results of 159 light 
curve solutions. Properties of contact binary stars were studied by using the catalogue data. 
As it is well known since Lucy's (1968a,b) and Mochnacki's (1981) works, primary components transfer their own energy to 
the secondary star via the common envelope around the two stars. This transfer was parameterized by a transfer parameter (ratio 
of the observed and intrinsic luminosities of the primary star). We proved that this transfer parameter is a simple function of 



' the mass and luminosity ratio. This newly found relation is valid for all systems except H type systems which have a different 

, relation. 

QO ' We introduced a new type of contact binary stars: H subtype systems which have a large mass ratio (q > 0.72). These systems 

I show highly different behaviour on the luminosity ratio - transfer parameter diagram from other systems and according to our 

. results the energy transfer rate is less efficient in them than in other type of contact binary stars. We also show that different 

' types of contact binaries have well defined locations on the mass ratio - luminosity ratio diagram. All contact binary systems do 

not follow Lucy's relation (L2/L1 = (M2/M1)'''-). No strict mass ratio - luminosity ratio relation of contact binary stars exists. 

o 

d 1. Introduction is convective. Although it explains the light curve shape and 
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the temperature equalization, however, contact binary star evo- 
Contact binary stars (or W UMa-type stars) consist of two j^^^^^^ -^^^^^^^ ^^^^^^^^^ ^^^^^^ photosphere re- 



dwarf stars whose spectral types are of F, G or K (only a ^^^^^ ^^^^ j^^^^^ ^^^^^^^ ^^^^^^^^^ ^^^^^^ ^^^^ discussed 
H few examples are known from earUer spectral types, and no M ^^^^^^ ^gge) and a recent review on them can be found 
^ spectral type contact binary star is known). Binnendijk (1965) Webbink (2003). Kahler (1989) summarized other possible 
pointed out that the components in a contact binary system have ^^^^^-^^ ^^^^^^^ ^-^^^^ ^^^^^ ^^^^^ ^^^^ ^^^^ ^p. 

nearly equal surface temperatures and luminosities in spite of p^.^^^^ observational facts. Recent detailed computations on 
their often highly different masses. If they are really dwarf energy transfer and internal structure were published by Kahler 
stars as suggested from their spectra, what is the mechanism (2002a 2002b 2003) 
which equalizes their temperatures and luminosities? The an- 
swer was given in Lucy's papers. The light curve characteris- 
tics of these binaries were succesfully interpreted by the con- 
tact model (Lucy 1968a,b), which simultaneously explains the 
shape of the light curve and the equal temperatures and lu- 
minosities of the components, hence this model is accepted 
for describing theoretically W UMa-type stars. The model as- 
sumes that both stars fill their Roche-lobe and therefore they B^^ed on the catalogue data, we investigated the efficiency 
are touching each other. Mass and luminosity is transfeiTed of *e energy transfer from the primary to the secondary star 
from the primary star to the secondary star through the nar- the mass ratio - luminosity ratio relation. These are outlined 
row neck between the components. In the model two main se- Sect. 3. 

quence stars are embedded in a common photosphere which According to the contact model, the energy generated in 
the cores of the components are redistributed in the common 



In order to develop our empirical knowledge about these 
stars, a catalogue on the results of the light curve (LC) solu- 
tions of contact binary stars - based on LC-solutions published 
formerly in the literature - was compiled and it contains the 
solutions of light curve of 159 systems. This catalogue is pre- 
sented in Sect. 2. 
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have another dependence on the mass ratio which would be the 
case if one looks at two main sequence stars in detached con- 
figuration (where Lz/^i = (M2/M1 )'*•<'). Lucy (1968a) found 
that the observable luminosity ratio is proportional to the ratio 
of the stellar surfaces (Lj/^i = (Mz/Mi)" ''^). We will show 
that the situation is more complex. 

Mochnacki (1981) assumed that the energy transfer rate 
from the primary star to the secondary one depends only on 
the mass ratio. Liu & Yang (2000) calculated the energy trans- 
fer rate and found its dependence on the mass ratio and on 
the evolutionary factor (which is defined as the ratio of the 
present radius and the zero-age radius of the primary compo- 
nent). Recently Kahler (2002a,b) examined this question with 
detailed numerical computations and he found the rate of the 
transferred luminosity to be variable in time. The transferred 
luminosity can vary within wide limits (see Fig. 2 of Kahler 
2002b) for a contact binary. (It should be emphasized that re- 
cent theoretical internal structure models of contact binary sys- 
tems are inconsistent, as was described in Kahler 2002b.) 

Kalimeris & Rovithis-Livaniou (2001) found that the ob- 
served rate of energy transfer is a function of the secondary's 
luminosity. We also examined empirically the energy transfer 
rate and found a simple relation between the luminosity and 
mass ratio and the amount of the transferred luminosity. 

As the referee of this paper pointed out the width of the 
neck is determined by the fill-out factor (which measures the 
degree of contact) and W type systems have thinner neck than 
A-type ones (see e.g. Mochnacki 1981; of course, our sample 
confirms this establishment) and have higher transfer rate and 
luminosity ratio than A-type systems (Figures 1 and 2). One 
can mind that the thickness of the neck determines the amount 
of transferred luminosity because through a thicker neck more 
luminosity can be transported, but this does not reaUzed. 

The reason of this paradoxon is not known yet. As a trend 
it is correct that the thinner the neck, the larger the rate of 
transferred luminosity which is exactly the opposite case than 
we expect. Note that Kahler's (2002b) numerical simulation 
showed that in general the degree of contact varies in phase 
with the transferred luminosity (compare Figures 1 and 2 in 
Kahler 2002b). The solution of this discrepancy between the- 
ory and observations requires further studies. 

2. The catalogue 

There are two recent catalogues on contact binary stars 
(Maceroni & van't Veer, 1996; Pribulla et al., 2003). The cata- 
logue of Maceroni & van't Veer (1996) lists 78 systems, while 
catalogue of Pribulla et al. (2003) contains 361 field contact 
binaries. The contents of these and our catalogues are demon- 
strated in Table 1 for comparison. 

If recent ephemeris was not available we repeated the 
ephemeris given in GCVS (Kholopov et al., 1998). Note that 
dimensionless surface potentials are generally assumed to be 
equal for the two components - the exceptions are noted in 
the catalogue. Gravity darkening exponents, albedos and limb- 
darkening coeflicients were generally fixed by the modellers 
with some exceptions, and these exceptions are noted in the 



Table 1. Comparison of contents of different catalogues of con- 
tact binaries. MV96: Maceroni & van't Veer 1996, PKT03: 
Pribulla etal., 2003. 



Content MV96 PKT03 This study 

photometric mass ratio XXX 
spectroscopic mass ratio XXX 
average fractional radii X X 

temperatures of the components XXX 
semiamplitude of V„rf of the primary star X 
magnitude of the O'Connell effect X X 

absolute dimensions (Ri,2,M\ 2, ^1.2 X 
of the components) 

system's angular momentum X 
ephemeris (epoch and period) X X 

number of available minima observations X 
information about period change X 
code which was used for light curve solution X 
inclination X X 

tillout factor X X 

spectral type and distance X 
range of variability X 
(nil + iut) sin'' ; X 
fractional luminosity of the primary X 
fractional luminosity of the third body X 
dimensionless surface potentials X 
spot parameters X 
gravity darkening exponents X 
albedos and limb darkening coefficients X 



catalogue. The catalogue together with its references is pub- 
lished in Table 2. 

For homogeneity, we primarily collected results of LC- 
solutions carried out by any version of the Wilson - Devinney 
code, however, in order to increase the sample results of mod- 
eling made by the BYNSIN Code (Vinko et al., 1992) was also 
included, but these are listed in a different table (see Table 3'). 

3. Relation between astrophysical quantities 

3.1. Subtypes of contact binaries 

In 1965, contact binary stars were divided into two sub-types: 
A-type systems (the larger star is the hotter one) and W-type 
systems (the smaller star is the hotter one) (Binnendijk 1965). 
Later Lucy & Wilson (1979) introduced the terminus of B-type 
systems which are systems in geometrical contact, but not in 
thermal contact and therefore there are high surface tempera- 
ture differences between the components. In this study we call 
B-type systems which has lOOOK or larger surface temperature 
difference between the components. Note that B-type systems 
are sometimes mentioned as PTC (Poor Thermal Contact) sys- 
tems (e. g., in Rucinski & Duerbeck 1997). 

' Tables 2 and 3 are available only electronically at the homepage 
of the Konkoly Observatory via the URL: ^http://www.konkoly.hu_ 
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Fig. 1. Mass ratio - luminosity ratio diagram of contact binary 
stars. Open squares: A, filled squares: B, filled circles: H/A (A 
subtype systems with q > 0.72), open circles: H/W (W subtype 
systems with q > 0.72). filled triangles: W subtype systems. 
Solid line is the /I = q^-^, the main sequence mass-luminosity 
relation. Dashed line is Lucy's relation A = q^-^. 
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Fig. 2. Transfer parameter vs bolometric luminosity ratio A. 
Symbols are the same as in Figure 1 . The solid line describes 
the minimum rate of transfer parameter. Note that most of the 
systems are on the this line. Systems which have q > 0.72 (H 
subtype systems) are far from this envelope. Dashed lines are 
due to different values of a (see text for further explanation) 
and the numbers show the corresponding values. 



We do not use the terminus 'E(arly)-type systems' which 
would mean the contact binaries of O, B and A spectral types. 
In Figures 1-4 they do not have any special position in con- 
trast to H-type systems. This conclusion is in agreement with 
Kahler's (1989) remark: from observational viewpoint there is 
no difference between early and late spectral type systems. 

In our sample we found 45 A, 13 B, 24 H and 77 W subtype 
systems. 
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Fig. 3. Transfer parameter jS corrected with the mass ratio vs 
bolometric luminosity ratio A. Symbols are the same as in 
Figure 1. 
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Fig. 4. Transfer parameter excess + A) ') vs mass ratio. 

Symbols are the same as in Figure 1. Solid line is the excess= 
0.52^4- 1 fit. 



3.2. Mass ratio - luminosity ratio relation 

Several systems were excluded from the next analysis because 
some authors did not publish all parameters of the light curve 
fitting. If fill-out factor was not given, we computed it with 
the BinMaker 2.0 software (Bradstreet 1992) from the mass 
ratio and surface potential, but other missing results could not 
be reproduced by us. That is why these systems were omitted 
from the sample. 

The temperatures of the components were known from 
modelings, therefore we could calculate the bolometric lumi- 
nosity ratios from the measured ratios in V band via A = 
(L2/Li)boi = (L2/Li)vl0°-4 <^'^i-^'^^'. Bolometric luminosity ra- 
tio vs mass ratio (q - M2/M1)) is shown in Fig. 1. (Bolometric 
correction was calculated from Flower's (1996) tables.) 

The subtypes of W UMa-stars are located at different re- 
gions in this diagram. W subtype systems have larger lumi- 
nosity ratio than A subtype ones at a given mass ratio. This 
is natural because the radius ratio is proportional to q^ '^^ at a 
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given mass ratio^, the temperature ratio TijTi is larger in W 
systems and hence the luminosity ratio is higher in W than in 
A systems. 

Both A and W subtype systems have larger luminosity ra- 
tios than any B subtype systems. A line represents the mass- 
luminosity relation if the components were main-sequence 
stars (in this case A - q^'^). Note that all systems are above this 
Une (with three exceptions: W Crv and LP UMa, and the B- 
type system HW Per). B subtype systems are located between 
A/W sytems and the A = q^-^ hne. 

Another interesting fact is that A-type systems are rela- 
tively rare objects at high mass ratios, but W-type systems show 
a completely opposite behaviour: they populate the region of 
higher mass ratios, and there are only few W-type systems be- 
low q - 0.3. This effect was discovered by Maceroni et al. 
(1985) while studying properties of W UMa stars based on a 
sample containing 42 systems. We confirm their result using a 
much larger sample. 

3.3. Energy transfer 

It is clear from Figure I that different mass ratio - luminosity 
ratio relations exist in the case of different subtypes. Systems 
also show remarkable diversity on the A-q plane. This is not due 
to the scatter of data. The luminosity ratio can be determined 
with an accuracy of 1-2% or better, while precision of mass ra- 
tio measured photometrically is generally better than a few per 
cent. In general, spectroscopicaUy and photometrically deter- 
mined mass ratios show a good agreement (Maceroni & van't 
Veer 1996, PribuUa et al. 2003). The diversity can be due to the 
dififerent rate of luminosity transferred from the primary to the 
secondary star. (Kahler, 2002b suggested such an efl'ect: in the 
same contact binary the transferred luminosity varies in time. 
At a given mass ratio this can cause diversity in luminosity.) 

We studied the energy transfer rate by the introduction of 
the transfer parameter. It was defined as 



minimum rate of the transfer parameter at a given luminosity 
ratio, so it is 



'^l, observed 
^1,ZAMS 



It is easy to show that p can be computed as 



or 



1 +q 



,4.6 



1+^0-92(^)4 



l+a^5.oi 



l+A 



(1) 



(2) 



(3) 



where a = (ri/r2)2" '". Transfer parameter was plotted against 
luminosity ratio (Fig. 2), and a good correlation was found with 
exception of all systems which have q > 0.72 (they are marked 
by a dififerent symbol in the Figures). The envelope is due to the 



^ This relation was found by Kuiper (1941), and it is a natural 
consequence of the Roche-geometry assumption in light curve mod- 
els. We also checked this relation applying the catalogue data and 
(R2/R1) = gO^f^o-ooa was found. 



^envelope — 



1 



l+A 



(4) 



It is interesting that systems with q > 0.72 are far from this 
envelope (hereafter we call these systems H-systems denoting 
high mass ratio systems) but other subtypes are close to it. 

To quantify the deviation from Eq. (4) we calculated the 
transfer parameter excess which was defined as the difference 
between j8 and the envelope given by Eq. (4). The excess was 
found to be a function of the mass ratio (Fig. 4) and a fit yielded 
that the he excess is proportional to 0.52(±0.02)^** '(+0.1). 

The systems should increase the luminosity transfer from 
the primary to the secondary in order to equalize the surface 
temperatures if the mass ratio is lower, but the situation is 
more complex than this simple picture. Note that Kalimeris & 
Rovithis-Livaniou (2001) found that the transferred luminos- 
ity is the function of the secondary's luminosity. Our results 
showed that the rate of the transferred luminosity is not related 
only to the luminosity ratio but it is the function of the mass 
ratio, too. 

In Figure 2 the yS-i curves for different a values (a = 0.25, 
0.5, 1.0, 2.0, 4.0 and 8.0) are also shown. (Note that a depends 
strongly on the surface temperature ratio: T1/T2 =1.11 corre- 
sponds to a « 8.) The a values of A, B and W subtype systems 
generally are close to 1 while H subtype systems have larger 
temperature ratio. If the mass ratio is close to 1, we would wait 
that the two components have very similar features and hence 
their surface temperature ratio (and a) is close to 1. In real- 
ity we found the opposite case: at large mass ratios the surface 
temperatures can be very different. 

From the definition of the transfer parameter it is clear that 
the amount of the transferred luminosity is AL = (1 - fj)Li. 
Substituting y8 = (1 +/1)"' -0.52q'^-^ we found that AL = (j^ - 
0.52q'^'^)Li. 

Figure 3 shows the corrected /3 (yScon = y6-0.52^'* ') against 
bolometric luminosity ratio. The correlation between them is 
very good confirming our conclusion that transferred luminos- 
ity is a function of mass and luminosity ratios. 



4. Summary 

The results of this research note can be summarized as follows. 

1 We compiled the catalogue of light curve solutions of con- 
tact binary stars. The catalogue contains LC-solutions of 
159 systems. 

2 We found that there is no strict mass ratio - luminosity ratio 
relation for contact binary stars. Such a relation was sug- 
gested by Lucy (1968a), but Kahler's (2002b) results in- 
dicated large luminosity ratio variations with small - prac- 
tically unobservable by light curve modeling - mass ratio 
variations. Kahler's model does not contradict our results. 

3 The energy transfer from the primary star to the secondary 
star was found to be depending on the mass ratio and the 
luminosity ratio. In H systems the energy transfer rate is 
less efficient than in other type ones at a given luminosity 
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ratio. We determined the amount of the transferred lumi- 
nosity from the primary to the secondary star and it was 
found a function of the mass ratio and the luminosity ratio. 
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